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We simplify and extend our previous model for the masses and mixing matrix of three Majorana 
neutrinos based on permutation symmetry 5*3 and the perturbations which violate this symmetry. 
The perturbations are arranged such that we get the smaller solar neutrinos mass difference at 
second order. We work out the corrections to the tribimaximal mixing matrix with the non-zero 
value for sis and the conventional CP-violating phase. It is shown that the results of the model are 
comparable with the global analysis of neutrino oscillation data. 

PACS numbers: 14.60.Pq, 12.15.F, 13.10.+q 



I. INTRODUCTION 

In the present paper, a certain approach to understanding neutrino mixing, proposed in [1| and further studied in 
0, 01, and [1| will be examined in more detail. This approach is based on assuming the permutation group S3 
symmetry for the lepton sector to be a starting point. The three neutrinos are arranged to be degenerate in this limit. 
Analogously to the SU(3) flavor and isospin breaking in QCD we add two breaking terms in different £3 "directions" 
with different strengths. Of course there have been many related treatments of this topic in the literature Q - [2{|. 

There are two new features compared to our earlier papers. First the order of the two perturbations will be arranged 
to correspond to the observed hierarchy between the "atmospheric" neutrino mixing and the "solar" neutrino mixing. 
Secondly the effect of "Dirac type" CP-violation will be included in addition to the previous Majorana type CP- 
violation. 
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II. BRIEF REVIEW AND NEW CHANGES 

In [l[ the "unperturbed" neutrino mass matrix was taken to be the S3 symmetric form, 
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al + /3d. 



(1) 
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a and j3 are, in general, complex numbers for the case of Majorana neutrinos while d is usually called the "democratic" 
matrix. 

This is diagonalized by, 



R T (al + f3d)R = 



a 
a + 3/3 
0a 



(2) 



where R has the so called "tribimaximal" form, 





R = ~ 75 — 75 — 7S 



(3) 
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FIG. 1: Isosceles triangle with angle ip between the equal length 2-vectors a and a + 3/3. 



To insure the initial neutrino degeneracy, the complex numbers a and a + 3/3 are arranged as shown in Fig. ([IJ . 
We introduce the notation, 



-i\a\e 



ii/>/2 



(4) 



The physical phase ip lies in the range: 
and | a | is related to /3 by, 



o < tp < it, (5) 
3/3 



1 1 2sm(V>/2)' 
In Q the effect of a perturbation matrix of the form, 



(6) 





A = | t u | , (7) 

x u t 

was considered with t and u real. 

In the present work we will first reduce the number of parameters by setting u = —t. This has the effect of eliminating 
the smaller "solar neutrino" mass splitting at the first order of perturbation. That, as well as CP-violation effects 
will be introduced at the next order of perturbation via the matrix, 



A' = ve'^ 




(8) 



Here v is real and & is a CP- violating phase. This represents a variation of the similar matrix A' used in the doubly 
perturbed model of [4] . 
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Hence the present model has only the three perturbation parameters t, v and <f>. t will be considered to be of first 
order strength while v will be taken to be second order. On the other hand the CP phase will not necessarily be taken 
small. The ve~ 1 ^ term will introduce both non-zero solar neutrino masses and Dirac type CP-violation. The main 
work will be to find the eigenvalues and eigenvectors of, 

M u = al + /3d + A + A'. (9) 



III. DIAGONALIZING THE NEUTRINO MASS MATRIX 



Since we are working in a basis where the zeroth order piece is diagonalized by the tribimaximal matrix R, we must 
bring to diagonal form the matrix: 



R 1 M V R 



( -i\a\e^/ 2 + \ve- 



V 



y5„ 





-i\a\e- 1 ^' 2 +3/3 




va ve -i<i> 



-i\a\e-^l 2 + 2t+\ve 



(10) 



Denoting the non-trivial 2x2, (13) sub-block by X, 



X = 



-i\a\e 



"A/2 i 3 



ve 



\(E ve -i<t> 



\ a \ e -i*t>/2 + 2t+ i ve -i<t> 



we diagonalize the hermitian quantity X^X with a unitary matrix U in the usual way, 

U^X^XU = Real, diagonal. 
This implies, as actually required to bring the Majorana type neutrino mass matrix to diagonal form, 

U T XU = diagonal, 



Explicitly, 



where to second order, 



U 



€ = 



\/3vsin(^ — <f>) 
4tsin($) 



1 e 

-e* 1 



t 



t 



\a\sin(^) \a\sin(-. 



(11) 

(12) 
(13) 

(14) 
(15) 



Here we used the approximations that t is first order, v is second order and v/t is, for consistency, first order. We 
trivially upgrade this to the 3x3 level as, 



satisfying, 



with, 



1 e 
Z = | 10 

-e* 1 



Z 1 R M V RZ = diag(d 1 ,d 2 ,d 3 ), 



di = -i\a\e-^/ 2 + -ve-^, 

d 2 = -i|a|e-^ /2 +3/3, 

d 3 = -i\ a \e~^ /2 +2t+ \e 



(16) 



(17) 



(18) 
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The neutrino masses uii are the magnitudes of the 



mi 

m 2 



| a- 1 
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3V ■ A\ 



lb v ,ib 
m 3 « \a\-2tsin(—)--sin(—-<f>). 

These mass parameters were made real, positive by the introduction of the phase matrix: 



P 



e~ lT 
e- u 




where, 
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The entire approximate diagonalization may be presented as, 

K T (al + f3d + A + A')K = 

where, 



mi 
m 2 
m 3 



is the full neutrino mixing matrix. Explicitly, 



K 



K = RZP 



-2 1 -2e 

V6 n/3 V6 

/F AT AT AT ' AT \ 1 , 



V6 i/2 V3 72 \/6 
\ * + Jl1 _L =1 + -L- 
\ \/6 V2 x/3 \/2 \/6 



(19) 



(20) 



(21) 



(22) 



(23) 



(24) 



where e is given in Eq. ([15 



IV. PARAMETER EVALUATION FROM EXPERIMENT 



First we have important information from neutrino oscillation experiments [30j-|37|. It is known that [38( 

A = m\- m? = (7.50 ± 0.20) x 10~ 5 eV 2 , 

B = \m\ - ml| = (2.32±g;^) x 10~ 3 eV 2 . (25) 
Also, constraints on cosmological structure formation yield [39l | a rough bound, 

mi + m 2 + m 3 < 0.3eV. (26) 

The two allowed spectrum types are: 

Typel (normal hierarchy) : m 3 > to 2 > mi, 

Type2 (inverted hierarchy) : m2 > rti\ > m 3 . (27) 
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From Eqs. ([19]) we write: 



m 2 — m\ = 3\a\vsin(— — (b), 

m 3 ~ m 2 — — 4|a|isin(-^) — |a|wszn(^- 



(28) 



For simplicity we will not try to determine all the parameters together. Rather we bound |a| at zeroth order, evaluate 
t\a\ at first order and evaluate v\a\ at second order. 

A. Zeroth order 

At this order the three neutrino masses are degenerate. With v = t = Eqs. (fHj|) and (|26p imply the bound, 

M < O.leV. (29) 

B. First order 

At this order the degeneracy between mi and is still preserved. With v = 0, t =^ one has, 

m 2 ~ m i = 



-4\a\tsin(^), 



(30) 



which implies, 



ta 



(31) 



This yields the following chart of t\a\ for several values of the Majorana-type phase ip and choice of the neutrino mass 
"hierarchy" . 







for normal hierarchy (eV 2 ) 


t\a\ for inverted hierarchy (eV 2 ) 


7T 


-5.8 x 1(T 4 


5.8 x 1CT 4 


tt/2 


-8.2 x 1(T 4 


8.2 x 10" 4 


tt/4 


-15.2 x 10 -4 


15.2 x 1CT 4 



TABLE I: t\a\ for different values of ib. 



C. Second order 

At second order the degeneracy between mi and m 2 is broken and v,t ^ 0. Using Eqs. ([28)) this determines v\a\, 

2 2 

mi — mi 

u|a| = \ ■ % \v (32) 

dsin{w — 0) 

We can plot v\a\ vs for various values of 'ib. As an example, the plot for ib = n is given in Fig. 2. 
The plots for other values of ib are similar. A more detailed comparison is furnished in Table II. 
Note that the ratio v/t of second order to first order parameters depends somewhat on the CP phase, <b but is less 
than 1/10 in magnitude. 



v|o| 
0.0002 r 




FIG. 2: Plot of v\a\ versus cj> for the choice ip = tt 





v\a\ for (j) — 


v\a\ for <f> = 7r/3 


v\a\ for (f> = 27r/3 


v\a\ for (j> = 7r 


v\a\ for <f> — 4n/3 


v\a\ for <f> — 5n/3 


TV 


2.5 x 10~ 5 


5.0 x 10" 5 


- 5.0 x 10" 5 


- 2.5 x 10" 5 


- 5.0 x 10~ 5 


5.0 x 10~ 5 


tt/2 


3.5 x 1(T 5 


-9.6 x 10~ 5 


- 2.6 x 10~ 5 


- 3.5 x 10" 5 


9.6 x 10~ 5 


2.6 x 10~ 5 


tt/4 


6.5 x 10~ 5 


- 4.1 x 10~ 5 


- 2.5 x 10" 5 


- 6.5 x 10~ 5 


4.1 x 10" 5 


2.5 x 10" 5 



TABLE II: v\a\ (in eV 2 ) for different values of ip and 

D. The parameter e and S13 

The complex mixing parameter e in Eq. (|15j) can be written as, 

-iS 



ee 



where, 



\/3vsin(^ — (j>) 
4f«n(f) 

V / 3wsm( : j — (j)) 
Usin(^) 



t 



tcos<j> 



The phase 8 is given as, 
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\a\sin(^) \a\sin(^ 

tcos(^)sincf) 
\a\sin{%)sin(% — 4 



tsin<fi 



\oc\sin(i — 4>) 



As mentioned above, we are keeping terms linear in t/|a|, v/t and u/|a|. Using this approximation we get, 

tsincj) 



sinS 



\a\sin(t — (j)) 



cosS ~ 1. 



Writing down elements of the matrix in Eq. (|24[) in terms of their magnitude and phases yields, 
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v sirup 
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(38) 



We would like to identify the transformation matrix K with a conventional form of the type: 



K eX pP 



— S12C23 
S12S23 " 



C12C13 

Cl2Sl3S23e Z 



Cl2Sl3C23e 



37 



S12C13 
C12C23 - Sl2Sl3S23e* 
-C12S23 - Si 2 Sl3C23e 



I'y 




(39) 



See, for example, Eq (25) of [J]. 

If we look at the matrix in Eq. (f3"T|). it seems like there are some extra phases, for example, the elements (23) and (33) 
have non-zero phases while (22) and (32) are without phases. In order to compare it with the above parametrization, 
we should take out these phases to bring it to the desired form. Let us multiply the matrix K in Eq. (|37[) by the 
diagonal matrices of phases L and Q from left and right. The resultant matrix will be identified with K exp , 



LKQ = K, 



exp • 



(40) 



LKQ 
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\/6 
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V3 
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V2j 6 
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V3 



1 
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iS 
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Ve 
M / 

76 / 



P' 



(41) 



Here. 



L = 



(- 


V 



1 








(42) 



Now the phase matrix L gets absorbed due to arbitrary phases for the charged leptons which will "sit" on the left 
side. What remains is the standard form when we identify P' = PQ as the Majorana phases. This enables one to 
easily compare with the three conventional CP conserving angles and the conventional Dirac phase. 

Comparing the matrices in Eqs. ([3TR) and (|4"Tj) . we first note that e in Eq. (|24l) is related to S13, 



Sl3 



-2|e| 
V6 ' 



(43) 



Also S23 and S12 can be read off: 



S23 



S12 



1 

71 



-0.58. 



V2 V6' 

Expanding S23 around its "tribimaximal value" as S23 = [s23]tbm + ^S23, one gets: 



(44) 



(S 23 ) 2 » ^ + ^2 ,5S 23, 

where SS23 measures the deviation of S23 from its tribimaximal value. 
From Eqs. (|43|) and p4|) . we find the interesting relation, 

s%3 = -2<5s 2 3- 



(45) 



(46) 



We should point out here that some similar results were obtained earlier in Q without the CP phase and with 
additional parameters. 
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Using the results of the global 3v oscillation analysis [40|, the numerical values of |e| and S13 are illustrated in Table 
IIIII for both normal (left table) and inverted (right table) hierarchies. It is encouraging that these values for |si3 1 are 
very well in agreement with the analysis of neutrino oscillation data. The best fit value for \s\%\ is 0.157 for normal 
hierarchy and 0.158 for inverted hierarchy which are also close to the predicted values at la, 2a and 3a levels. 

Note also that the magnitude of parameter e and S13 are only determined from the deviation in the mixing angle 
S23 so the corrections to S23 are fixed as our input. The mixing angle S12 does not get any correction in this model 
but it is close to the best fit value of 0.56 from the data. 





la 


2a 


3a 


la 


2a 


3a 


{\5s 2 'i\)data 


0.061 


0.10 


0.11 


(\8s 2 3\)data 


0.069 


0.09 


0.11 


N 


0.149 


0.24 


0.27 


N 


0.169 


0.22 


0.27 


1*13 1 


0.122 


0.20 


0.22 


l*13| 


0.138 


0.18 


0.22 


(Sis)data 


0.157 


0.16 


0.16 


(Sl3)data 


0.158 


0.16 


0.16 



TABLE III: Numerical values of |e| and S13 for normal (left) and inverted (right) hierarchies. 

Finally the phase 6 in Eq. (j3"5)l can be estimated using Table |U Some values are given in Table IIVI for both mass 
hierarchies. The best fit value for the phase S in the global analysis of neutrinos is 0.80-7T for normal hierarchy and 
— 0.03-7T for inverted hierarchy. It can be seen that our values are close to the values for inverted hierarchy but an 
order smaller for normal hierarchy. 





5 for <j) = 


5 for <j> — 7r/3 


5 for <f> = 2tt/3 


S for (f> — ty 


S for <j) = 4?r/3 


5 for cj> = 57r/3 


7T 





0.037T 


- 0.037T 





0.037T 


- 0.037T 


tt/2 





- 0.097T 


- 0.027T 





- 0.097T 


- 0.027T 


7r/4 





- 0.077T 


- 0.047T 





- 0.077T 


- 0.047T 



TABLE IV: S for different values of ip and (j> for normal hierarchy (opposite sign for inverse hierarchy). 



V. SUMMARY AND DISCUSSION 

In this work, we studied the model [l[ , [|| and Q for neutrino masses and mixing based on the permutation symmetry 
S3 in more detailed and simplified way. In Sec. II, we presented the outlines and differences from previously studied 
models. Our starting point was the full S3 invariance, with parameter |a|, where we imposed the unperturbed neutrino 
mass degeneracy by introducing a Majorana type phase ip. Next, we added a perturbation symmetric under the (23) 
subgroup of S3. We chose the parameter t for this perturbation in a way to eliminate the smaller "solar neutrino" 
mass difference. The second perturbation with complex parameter ve~ 1 ^ , living in the (12) subspace, was added in a 
similar way with (f> related to Dirac type phase. 

In Sec III and IV, we carried out the diagonalization of the resultant complex symmetric matrix and determined 
the parameters of the model. We used the approximation in which we kept terms of the type i/|o;|, v/t and ?V|a|. 
The ratio v/t is effectively of first order and is less than 1/10 in magnitude. The solar neutrino mass difference came 
out to be proportional to the second order parameter v and the atmospheric neutrino mass difference to the first order 
t as well as second order v. Both mass differences were also dependent on the Majorana phase ip and the Dirac phase 
4>. For the mixing angles, the corrections had a piece proportional to the complex parameter e which is effectively of 
the first order parameter v/t. 

Important motivations for this model were the non-zero value of S13 and the Dirac CP phase. We found that 
the mixing angle | s 13 1 was related to the deviation of S23 from its tribimaximal value which we used as one of our 
experimental inputs. The predicted values for | S13 1 were close to the numbers from data (see Table ILTT)) . Our results for 
the CP-violating phase S were well in agreement with the suggested value from global analysis of neutrino oscillation 
data 40] for inverse hierarchy and an order less for normal hierarchy (see Table IIV[) . The mixing angle S12 was 
unchanged from its tribimaximal value by the perturbation in this model. 

The lcptonic factor describing the amplitude for the neutrinoless double beta decay of a nucleus: (A, Z) — > (A, Z + 
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2) + 2e comes out to be, 

\m ee \ w -|a| a/5 + 4cosV 

It may be noted that |m ee | depends also on the Dirac CP phase <fi but this effect is negligible since it appears at 
second order. Some predictions for \m ee \ at zeroth order for typical values of ip are listed in the following chart. 

ip: n, 7r/2, 7t/4 

\m ee \ : 0.033 eV, 0.074 eV, 0.093 eV 



1 \a\ 



bcos^sinq 



5 + Acosip 



(47) 



The current experimental bound on \m ee \ is < 0.36 eV [41| . It can be seen that \m ee \ gets larger for smaller values 
of the Majorana phase, ip, for degenerate neutrino masses. 

We have not discussed the role of the Higgs sector in this paper but some technical details were discussed in [l[ . 
We have also assumed that the diagonalization of the charged lepton mass matrix makes a negligible contribution to 
the lepton mixing matrix. Of course, at a later stage, the contributions of the charged leptons should be included at 
the expense of additional parameters or assumptions. Further results on this model will be given elsewhere. 
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